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In this paper we show that electron-positron pairs can be pumped inexhaustibly with a constant 
production rate from the one-dimensional well potential with oscillating depth or width. Bound 
states embedded in the the Dirac sea can be pulled out and pushed to the positive continuum, 
become scattering states. Pauli block, which dominant the saturation of pair creation in the static 
super-critical well potential, can be broken by the ejection of electrons. We find that the width 
oscillating mode is more efficient that the depth oscillating mode. In the adiabatic limit, pair 
number as a function of upper boundary of the oscillating, will reveal the diving of the bound 
states. 
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I. INTRODUCTION 

Since Einstein’s relativistic theory tell that the mat¬ 
ter can convert into energy, the possibility of converting 
energy into matter, i.e., the electron-positron pair, as 
predicted by Dirac in quantum electrodynamics [ 1 ], has 
attracted a great deal of interest [2] . In presence of a static 
and uniform electric field, the quantum electrodynamic 
(QED) vacuum may break down and decay into electron- 
positron pairs due to a quantum tunneling effect [4-6]. 
The critical Schwinger field is = rn^e? j (|e| fi,), which 
can accelerate the electron to an energy of the order of 
its rest mass on its Compton wavelength Ac = h/mc, 
where m is the electron mass. Starting from the works of 
Brezin and Popov et.al [7-9], the Schwinger mechanism 
was generalized to time dependent fields [10-17] , where 
another mechanism may be responsible for the pair cre¬ 
ation. If the frequency of the alternating field exceeds the 
gap 2mc^, electrons in Dirac sea can transit to positive 
states and pairs are triggered. Experimentally, pairs can 
be generated by the relativistic heavy-ion collisions [18] or 
the collision of an intense laser pulse and a 46 Gev elec¬ 
tron beam[19[, but pairs created from pure laser light has 
not been observed until now. 

Recently, various numerical approach were developed 
to deal with the time dependent Dirac equation[20- 
22], the Klein paradox[23, 24], the Zitterbewegung[25[, 
and the pair production process [26-31]. The one¬ 
dimensional well potential, specially, for its simplicity, 
is studied extensively [27-29, 32, 33] . The super critical 
well potential has bound states embedded in the nega¬ 
tive continuum can cause spontaneous electron-positron 
pair creation. Theoretical investigations are expected 
to make the physics of the creation clear and predict 
a higher generation rate. However, The Pauli exclu¬ 
sion principle will block further creation once the bound 
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states are occupied, resulting a asymptotic saturation 
behavior[27, 28, 33]. Motivated by this requirement and 
a better understanding of the pair creation process in the 
one-dimensional well potential, we examine the pair cre¬ 
ation in a well with its width or depth oscillating. By 
oscillating the width and depth, the transfer channels for 
population are opened and closed alternately. The elec¬ 
trons confined in the well will be released and the Pauli 
block become invalid. This can lead to a non-vanishing 
production rate, which means that pairs can be pumped 
inexhaustibly form the well. 

This paper is organized as follows. In Sec.II we present 
the model and the numerical method we employed. The 
well potential is set to be oscillating in two modes, the 
width oscillating mode and the depth oscillating mode. 
The energy spectrum is plotted as a function of the width 
or the depth. In Sec.III. we discuss the pair production 
process in both two modes. The time evolution of pair 
number, spacial density and pumping rate are studied. 
We also investigate the adiabatic limit of the oscillating. 
In the last section we give a brief summary. 

II. MODEL AND METHOD 

A. Model: one-dimensional well potential with 
oscillating depth or width 

In one dimension, the time evolution of the Heisenberg 
field operator T (z, t) is given by the Dirac equation ( 
without spin, for simplicity in this paper ) [34] 

t) = • Pz + + V (z, t)] ^ (z, t). (1) 

(Ti (T 3 are Pauli matrices, c is the speed of light in vac¬ 
uum, V (z, t) is the external potential. The atomic units 
([a.u.]) is used in this paper: m = h = e = l, c=l/aKi 
137.0359991, a is fine-structure constant, Compton wave 
length of electron is Ac = 1/c. The Hamiltonian of the 
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Figure 1. The energy spectrum of the total Hamiltonian 
as a function of the width or the depth of the potential, 
(a), Vo = 2.53c^, as W increasing, the bound states dive 
into the Dirac sea at W — 2.79, 5.51,8.21...(in units of Ac, 
the electron Compton wavelength), (b), W = lOAc, as 
Vo increasing, the bound states dive into the Dirac sea at 
Vo = 2.05,2.19,2.38,2.62,2.87,3.15,3.43,3.73, ...(in units of 
c^). 


system is H = \ccti ■ + V (z, t)]. We define 

the potential as 
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D is the width of potential edge (a measure of the width 
of the electric field), and we set D = 0.3Ac- The numer¬ 
ical box size is set to L = 2.5. 

The potential width W (t) and the depth Vq (t) ( 
positive, but note that the potential V (z, t) is nega¬ 
tive in the center and zero elsewhere) are set to two 
modes: (1) the W-oscillating mode: Vq is constant, 
W{t) =Wi + ^ {W 2 - Wi) [1 + sm{u;w{t) - 7r/2)] ; (2) 
the V-oscillating mode; W is constant, Vq (t) = Vi + 
i (V 2 — hi) [1 + sin(a;y(t) — 7r/2)]. In this paper we as¬ 
sume Wi = 0 and Vi = 0, then W {t) ( or Vq {t)) varies as 
a sine function between zero and its upper boundary W 2 
( or V 2 ), and the turning on process is from zero point 
with first order derivative equal to zero. In the following 
numerical simulation, we choose the total evolution time 
to be the period ( Tw or Ty ) of the oscillating W (t) or 
V (t) multiples an integer, to make the potential turning 
off is finished with first order derivative equal to zero. 

The numerically energy spectrum of the total Hamil¬ 
tonian of finite-size (length, for one dimension) are pre¬ 
sented in Fig. 1 for varying W and Vq, which can be a 
schematic of the real one-dimensional system. In Fig. 1, 
we show the critical width or depth, and the behavior 
of ’diving’ of the bound states into the negative contin¬ 
uum. For example, if Vq = 2.53c^, there are bound states 
embedded and then pair can be spontaneously triggered 
only when W > 2.79Ac- 


B. Method: the numerical quantum field 
theoretical approach 

In recent years, numerical quantum field theoretical 
approach [31] has been established to overcome the sin¬ 
gle particle picture described by quantum mechanics and 
the mathematical difficulty of quantum electrodynamics. 
In this section we will briefly review this method and de¬ 
scribe that how do we deal with the model in this paper. 

The field operator can be expressed in terms of the elec¬ 
tron annihilation and positron creation operators as [31] 


p n 

= ^&p(t)Wp(z) + ^dl(t)W„(z), (4) 

p n 

in which Wp(„) (z) = (z|p(n)) is the solution of the filed- 
free Dirac Hamiltonian ( V (z, t) = 0), (z, t) = 

{z\p{n) (t)) is the time dependent solution of the Dirac 
equation ( 1 ), and the term X)p(n) denotes the summa¬ 
tion over all states with positive (negative) energy. The 
eigenstates of the filed-free Hamiltonian are 


Wp(z) 

Wniz) 


y/E + d^ 


y/VnE \_sign (p) y/E - c^_ 

-sign (n) y/—E — 


y/—4TrE 


y/-E + c2 


(5) 

( 6 ) 


where Ep = y/cd p^c^, and E^ = —y/c^ respec¬ 

tively. The time dependent single particle wave function 
Wp(„) (z, t) can be got by introducing the time-evolution 

operator U (t 2 , ii) = Texp dt'H (t')^, 


%(„)(z,t) = f7(t, t = 0)Wp(„)(z), (7) 


where T denotes the Dyson time ordering operator. In 
this paper, we use the numerical split operator technique 
[ 20 , 21 ], then 

W{t + dt) « (t) 

= e-*f -b O (dt^) , ( 8 ) 


with 


Ho = C(Ti ■ p^ + C^(T3, (9) 

H, = V (z, t). (10) 

Practically, since the derivation ( the momentum opera¬ 
tor ) can be implemented by replacing the operator 
with its value in momentum space, the evolution op¬ 
eration has the following form 






























III. PUMP ELECTRON-POSITRON PAIRS 
FROM THE WELL POTENTIAL 
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jF-l[cos {(j)) 


—i sin {(j)) 


ai- kz + cas 


+ kl 




—i sin (y (t) dt)]W (t), (12) 


W (t) = [cos (V (t) dt) 


where kz is momentum, </> = ^+ k^, and F 
is Fourier transformation (inverse Fourier transforma¬ 
tion) . 

Then, after the time dependent field operator 'h (z, t) 
can be calculated, the number, the spacial distribution of 
electrons created from the vacuum ( defined as bp jjuac) = 
0 , dn jjuac) = 0) are obtained from the positive part of 
the field operator. 


iV®*' (t) = (uacjj (a;, t) {x, t) jjuac) 

pn 


Nf (t) = ^ 

n 


{t)Wp{z) 


(13) 

(14) 


where Up„ (t) = {Wp (z) |1F„ {z, t)) = 

J dxWyz)Wn {z, t). The pair number N (t) is 
equal to the electron number A®*' (t). 

The spacial distribution of the created positrons can 
be written as 


Ny (t) = 5] 


Upn it) W„ (z) 




(15) 


The total positron number Np°- (t) is equal to the elec¬ 
tron number A®*' (t). 

We can also get it from the negative part of the field 
operator by compute the number and spacial distribu¬ 
tion of the holes. In this paper we use this expression ( 
Eq. (15)) to reduce the computational cost, because Upn 
has been calculated in Eq.(13). Furthermore, we can ne¬ 
glect the larger part of the momentum ( y/Wdy+l? is 
far greater than V and w) in the numerical simulation, 
for its contribution to the matrix element Upn (t) is very 
small. In the following, the number of spatial points is 
Nz = 2048, and we only take Np = 1024 discrete mo¬ 
mentum into account. 

Based on the projection of the field operator onto the 
the field-free electronic states in this method and the def¬ 
inition of electron and positron in the Dirac hole theory, 
in this paper we will present physical quantities for all 
time and focus on the moments when the field is absent. 


For a well potential of depth Vb, if Vq < 2c^, the posi¬ 
tive continuum and negative continuum can not overlap. 
But for a super-critical depth, Vq > 2c^, the domain 

— Vo < E < —c^ exist, and bound states in the well 
are possible (which we call ’ bound states embedded in 
the negative continuum ’ ) : their wave function do not 
decrease exponentially out the well, but join a continuum 
wave of the same energy E < —c^ out the well, hence the 
wave function has a non-zero probability outside. An 
empty bound state will spontaneously be occupied by an 
electron (two, if the spin is considered) from the filled 
Dirac sea, and the hole (identified as positron) will travel 
away from the well to infinity [34]. This is the picture 
of spontaneous creation of electron-positron pair. For 
a static well potential, electrons will fill the embedded 
bound states, and the Pauli principle will prevent fur¬ 
ther pair creation. The number of pair created should be 
the number of bound states which meet these conditions. 

For a time dependent potential, the situation is more 
complicated. In paper [28], the effect of open and close 
a pair-creation channel was studied. The well depth is 
fixed at Vb = 2.53c^, while the width W varies between 
Wi = 4.55Ac and W 2 = 6.15Ac. For W = 4.55Ac, 
there is one bound state embedded in the Dirac sea, and 
there are two for W = 6.15Ac. After enough time for 
saturation, the pair number will increase as one more 
channel is opened, but do not decrease as one of the two 
channels is closed. The reason is that the annihilation of 
the pair need the electron and positron to be in the same 
place, which is not satisfied because the electron remains 
in the well while the positrons have left the creation zone 
and escaped to the opposite direction. 

Naturally, one can propose that if the channel is 
opened and closed periodically, can this mechanism will 
lead to a continuously pair creation? Moreover, for fixed 
W and varying Vb, since the diving behavior is similar 
(Fig. 1), will something similar happen? Motivated by 
these questions, we construct two oscillating modes as 
described in Sec. IF A. Results and discussion are as 
follows. 


A. time evolution of pair number 

Using the method presented in Sec. II. B, we graph the 
time evolution of the pair number defined as Eq.(13) for 
both W-oscillating and V-oscillating mode in Fig. 2. The 
width frequency cuw and depth frequency ww are in units 
of c^, and their values are assumed to be relative low, 
comparing to the gap 2c^, so that the photon absorption 
mechanism is not valid. The total time is 1207r/c^ « 0.02 
and the period is Tiy = 2tt juiw or Ty = 2tt juiy. The 
dot line represent the time t ^ L/ (2c) « 0.009 when the 
particles arrive the boundary, z = ±L/2 = ±1.25. Since 
IVi = 0 and Vi = 0, if the time is an integer multiples of 
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Figure 2. The time evolution of the total number of pairs for 
both W-oscillating and V-oscillating mode, (a), W-oscillating 
mode, W 2 = 10Ac,Vb = 2.53c^; (b), V-oscillating mode, 
V 2 = 2.53c^, W = lOAc- The frequency cJw and cvv are 
in units of c^. The dash line represent the time t = 0.009 
when the positrons arrive the boundary, 2 : = ±LI2 = ±1.25. 
The triangles denote pair number when the field is absent. 
The dot line just link these triangles. 


the period ( Tw or Ty ), the system Hamiltonian degen¬ 
erate to a field free one. The triangles in Fig. 2 denote 
the pair number when the field is absent. 

W-oscillating mode: In Fig. 2 (a), we illuminate 
the total number of pairs as a function of time for wp/ = 
0.1/6c^, 0.2/3c^, 0.3c^, 0.6c^. The depth Vq is fixed at 
Vo = 2.53c^. The width W varies between Wi = 0 and 
W 2 = lOAc, corresponding zero and three bound states 
embedded. When W = W 2 , there are also eight bound 
states exist in the gap, which can be associate with the 
pair creation[27]. 

When uJw = 0.1/6c^, the width W can only finish one 
cycle in the total time 1207r/c^. JV begin to arise before 
t = 3.57X10“^ , corresponding W (t) = 2.79Ac, when the 
first bound state dive into the negative continuum. The 
reason is the non-adiabatic varying width, and JV will 
begin to arise precisely at the time when W (t) = 2.79A(y 
in the adiabatic case (ww 0 , see the discussion be¬ 
low) . JV increases as more bound states dive in, and reach 
its maximum N = 2.89 at t = 1.37 x 10“^ , between 
t = 1.28 X 10“^ and 1.47 x 10“^, at which time the third 
and the second bound state were pulled out the Dirac sea. 
Undergoing the particle-antiparticle annihilation, N de¬ 
creases but remains an appreciable value N = 2.85 at the 
end. In the latter half of this cycle, the embedded bound 
states depart from the Dirac sea, return to the positive 


continuum, and become scattering states. The released 
positrons are reflected by the numerical box boundary, 
come back to the interaction region and will affect the 
pair generation after. Though the effect is weak when 
ujw = 0.1/6c^, it is non-ignorable when, i.e., loy/ = 0.3c^ 
(see Fig. 3 for details). 

For wiv = 0.2/3c^ and ujy/ = 0.3c^ , W can finish four 
and eighteen cycles in the total time and the pair number 
are N = 6.49, N = 21.4 at the end. For t < 0.009, W can 
finish one and eight cycles, respectively. In each cycle, the 
positrons are repulsed by the electric field to the infinity 
once they were generated, while the electrons are limited 
in the well when the field is strong enough and extruded 
out as the well is turning off, avoiding the inevitable Pauli 
block in the non-varied static well construction. The non- 
synchronous ejection prevent the annihilation and lead to 
a high production rate. 

The next cycle starts from field free and is independent 
on the previous cycle. In Fig. 2, the dot line link the 
triangles which denote the pair number when the field is 
absent. We can find that the pair generation denoted by 
the dot line is linearly depend on time for low frequency 
uJw, for t < 0.009. If the system length L is infinite and 
there is no reflection at the boundary, the pairs can be 
pumped inexhaustibly with a constant production rate 
from the well. Even for ww — 0.6c^, there is nonlinear 
effect at the beginning, the generation rate become stable 
soon. 

Due to the finite period T\y and the bound states in 
the gap, particle generation and ejection process is not 
monotonic with the increase of the frequency UJw^ see 
Fig. 2 (a). However, ignoring the reflection, if the W- 
oscillating frequency uJw is very small, we can expect a 
linear dependent of final pair number on the frequency. 

V-oscillating mode: The number of pairs N as a 
function of time are presented in Fig. 2 (b), for ujy = 
0.1/6c^, 0.2/3c^, 0.3c^ and 0.6c^. The width W is fixed 
at VF = 10Ac, while the depth varies between Vi = 0 and 
V 2 = 2.53c^, corresponding zero and three bound states 
embedded. There are also eight bound states exist in the 
gap when Vq = V 2 - 

For uJw = 0.1/6c^, the first bound state dive in at 
t = 7.20 X 10“^, at which time there are already N = 
8.83 X 10“^ pair generated. The first bound state depart 
the negative continuum after the second and the third 
one, at t = 1.29 x 10“^, when N reach its maximum 
N = 1.81. Finally, there are N — 1.74 pairs survived at 
t = 1207r/c^. For uiy = 0.2/3c^, 0.3c^, 0.6c^, the pair 
number at the end are N = 2.21, 2.56, 3.78. 

Instead of pulling and pushing the walls of the well 
in W-oscillating mode, in this mode it is the rising and 
falling bottom of the well that control the bound states 
diving in and departing from the negative continuum. It 
is also the non-synchronous ejection of the positrons and 
electrons which dominant the pumping process. 

The dot line here indicate a linear relation between the 
pair number and time. The final number is not mono¬ 
tonic depending on the frequency cuy, and we can also 
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Figure 3. For W-oscillating mode, uJw = 0.3c^, the three 
dimensional diagrams for entire time and the waterfall figures 
for held free moments (the time indicated by triangles in Fig. 
2 (a)), for electron spacial density (a, c) and positron spacial 
density (b, d). The thicker curve in sub-hgure (c, d) mark the 
last cycle before the positron arrive the boundary. The well 
potential V{z) with Vo = 2.53c^, W = lOAc, are included 
on the bottom for comparison. All other parameters are the 
same as Fig. 2 (a). 


Figure 4. For V-oscillating mode, ujv = 0.3c^, the three di¬ 
mensional diagrams for entire time and the waterfall hgures 
for held free moments ( the time indicated by triangles in Fig. 
2 (b)), for electron spacial density (a, c) and positron spacial 
density (b, d). The thicker curve in sub-hgure (c, d) mark the 
last cycle before the positron arrive the boundary. The well 
potential V{z) with Vo = 2.53c^, W = lOAc, are included 
on the bottom for comparison. All other parameters are the 
same as Fig. 2 (b). 


expect a linear dependent of final pair number on ujv 
when u>v is very small. 

Note that although the two modes has the same begin¬ 
ning and ending parameters, the generation rate in the 
W-oscillating mode is much higher. 


B. time evolution of spacial density 

In last section we discussed the total pair number as 
a function of time in different oscillating frequency for 
both modes. To show the pumping process explicitly, we 
compute the time evolution of spacial density of electrons 
and positrons ( Eq. 14 and Eq.15) for uJw = 0.3c^ and 
ojv = 0.3c^ respectively. 

In Fig. 3, for W-oscillating mode, uw = 0.3c^, we 
plot the the time evolution of spacial density of electrons 
and positrons (sub-figure (a) and (b)). Specially, for the 
moments when the field are zero, denoted by the trian¬ 
gles in Fig. 2 (a), these quantities are plotted in the 
waterfall figures. Fig. 3(c, d). For V-oscillating mode, 
uJv = 0.3c^, similar diagram are presented in Fig. 4. 
For comparison, the well potential V{z) with wide and 
depth equal to the upper boundary of the two modes, 
Vq = V 2 = 2.53c^, W = W 2 = lOAc, are included on 
the bottom. These figures clearly show the process how 
are the particles pumped from the well and spread in the 
numerical box. 

Since cuw = 0.3c^, The period of the width oscillating 
is Tiy = 1.12 X 10“^. Before positrons arrive the bound¬ 
aries, the width can finish eight cycles. If we detect the 


particle population at the boundary, we can find that 
positrons arrive the boundary first, at t = 9.15 x 10“^, in 
conformity to the estimation Lj (2c) = 9.12 x 10“^. Elec¬ 
trons arrive the boundary at t = 1.02 x 10“^, about one 
period (Tw or Ty ) later than the positrons. We can see 
that the particles reflected by the boundary come back 
to the interaction region, and may cause non-ignorable 
effect, i.e., the non-linearity of the last three triangles in 
the dot line in Fig. 2(a), uJw = 0.3c^. 

Comparing with the rising and falling bottom of the 
well, more work is done by the wall of the well in the 
case of opening and closing the well. In the W-oscillating 
mode, the wavefront of the particles are more abrupter 
and regular. In energy space, higher energy modes are 
excited, and the spectrum show periodic structure with 
0.3c^ between each peak. In the V-oscillating mode, elec¬ 
trons are lifted and released naturally. Less work is done 
and only low momentum mode are excited, the rate of 
electrons in the well region (—5Ac < 2 < 5Ac) is more 
larger. Also, We can see the absent of interferences both 
in or out the well, as discussed in [27]. 


C. time evolution of pumping rate 

In the V-oscillating mode, it turns out that the elec¬ 
trons is more inclined to gather in the well region (de¬ 
fined as —5Ac < z < 5Ac) than it in W-oscillating 
mode . We can integrate the spacial density N (z) in 
this region and get the particle number in the well, 
^et(po.)(^) ^ j5Xc^ For the pumping pro- 
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Figure 5. For W-oscillating mode {ujw = 0.3c^, sub-figure 
a, c) and V-oscillating mode {u>v ~ 0.3c^, sub-figure b, d), 
particles in the well {Nin) and the pumping rate Nout/N as a 
function of time. The triangles denote the time when field Is 
absent and the dot line link them. The blue triangles denote 
electron and the red denote positron. All parameters are the 
same as Fig. 3 and Fig. 4, respectively. 




Figure 6. The final number of pairs created after one cycle 
as a function of the upper boundary of the oscillating width 
and depth, (a), W-oscillating mode, Vo = 2.53c^; (b), V- 
oscillating mode, W = lOAc. The total time T is chosen 
equal to the oscillating period. W 2 is in units of Ac, T is in 
units of 1/c^, and V 2 is in units of c^. 


cess in last section, are graphed in Fig. 5(a, 

b). In W-oscillating mode, as time increasing, in¬ 
crease to a constant 1.60 quickly, while Nf°' to a constant 
0.36. But in W-oscillating mode, keep increasing 
while Ad^'keep zero. The reason is positrons can be gen¬ 
erated in the well region in W-oscillating mode, while the 
wall (the electric field) shut the door upon positrons in 
V-oscillating mode. 

In a pumping process, the pumping rate is vitally im¬ 
portant and can be defined as a{t) = Nout/N, where 
Nout = N — Nin, as shown in Fig. 5(c, d). In both 
modes, at the end of the first cycle, when t = T\y or 
Tv, nearly all the electrons are limited in the well region, 
while positron are ejected. In V-oscillating mode, since 
all the generated positrons are ejected and kept out of 
the well, the pump rate directly become 1. For electron 
in the V-oscillating mode, or electron and positron in W- 
oscillating mode, in the long time limit, a(t) come to 1 
as 1 — (3/1, where (3 depends on the saturation number of 
particles in the well and the number of particles can be 
generated in each cycle. 


D. The adiabatic limit 

In Fig. 2, for uJw = 0.1/6c^ and coy = 0.1/6c^, there 
are N = 2.85 and N = 1.74 pairs survived at the end 
t = 1207r/c^. We have proposed that in low frequency 
limit, the pairs survived finally should equal to three, 
the maximum number of embedded bound states swept 
in one cycle of each mode. In Fig. 6, ignoring the reflec¬ 
tion, for each frequency, the total time is chosen equal 
to the oscillating period for both modes, so that the os¬ 
cillation can only finish one cycle. The final number of 


pairs survived as a function of the upper boundary of the 
oscillating width ( W 2 ) and depth ( V 2 ) are presented. 

In the adiabatic limit, a sub-critical well potential can 
not trigger pairs. As the width or depth increasing, the 
bound states in the gap dive into the negative continuum 
successively, the potential become super-critical. Pairs 
can be generated and saturated to the number of em¬ 
bedded bound states. However, as the width and depth 
decreasing, bound states depart the negative continuum 
successively and the generated pairs can not annihilate 
because of the non-synchronous ejection. Finally, the 
number of pairs survived in this cycle is equal to the 
maximum number of bound states embedded. This max¬ 
imum number is a function of the upper boundary of the 
two oscillating cycle ( W 2 or V 2 )• 

For a low frequency, the curve which indicate the fi¬ 
nal number of pairs vs. W 2 or V 2 is like a flight of 
stairs. As the frequency become lower, the rising edge 
of the stairs become more sharper. As shown in Fig. 
6 , in the limit wvVj 0 , the the rising edge of 

the stairs will precisely locate at the points where the 
bound states dive into the negative continuum. These 
points are W = 2.79, 5.51,8.21...(in units of Ac), and 
Vo = 2.05, 2.19, 2.38, 2.62, 2.87, 3.15,3.43,3.73, ...(in units 
of c^), as illuminated in Fig. 1. 

The gaps between bound states in the positive and neg¬ 
ative continuum in V-oscillating mode are smaller than 
that in W-oscillating mode. To achieve a quasi-adiabatic 
( finite Tw or Ty) simulation, Ty should be larger than 
Tyv to build a similar stairs. 

Now, if the two quasi-adiabatic oscillating cycle re- 
peate periodically, we can expect a linear increasing pair 
number, i.e., for Fig. 6(a), W 2 = 7Ac, the final pair 
number will be 2 times the number of the cycles. 
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IV. SUMMARY 

In this work, we have constructed a toy model, one¬ 
dimensional well potential with its width and depth oscil¬ 
lating, and studied the electron-positron creation. Since 
the bound states diving behavior in the energy spectrum 
are similar when sweeping the depth or width, the phys¬ 
ical process are similar in these two modes. We find that 
the non-synchronous ejection of particles prevent the par¬ 
ticle annihilation, break the Pauli block effect in a static 
super-critical well potential, and lead to a high constant 
production rate. The width oscillating mode can deliver 
more energy to particles and is more efficient in pumping 
pairs than the depth oscillating mode. The time evolu¬ 
tion of spacial density illustrate the particles pumping 
from the well and the spreading of them in the numerical 
box. In a quasi-adiabatic case, pair number as a function 
of upper boundary of the oscillating, will reveal the div¬ 
ing of the bound states. This can be expected to detect 
the energy structure of a complicated potential. 

In order to reduce the computing cost, we neglect the 
larger part of the discrete momentum in the numerical 


simulation. On the other hand, with the same number of 
discrete momentum, the number of spatial points can be 
larger to describe the details of the potential. The sim¬ 
ulation in this paper is done on a personal stand-alone 
computer. In this algorithm, the time evolution of each 
negative eigenstate can be done on a single CPU, hence 
the computation can be paralleled easily. Further more, 
if the second order spatial derivative in the Hamiltonian 
are done by finite difference approximations instead of 
Fourier transformation[22], more lager one-dimensional, 
even two dimensional system can be simulated through 
paralleling the algorithm on memory shared parallel com¬ 
puters. 
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